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ABSTRACT 

Microlensing events are usually selected among single-peaked non-repeating light 
curves in order to avoid confusion with variable stars. However, a microlensing event 
may exhibit a second microlensing brightening episode when the source or/and the 
lens is a binary system. A careful analysis of these repeating events provides an in- 
dependent way to study the statistics of wide binary stars and to detect extrasolar 
planets. Previous theoretical studies predicted that 0.5 - 2 % of events should repeat 
due to wide binary lenses. We present a systematic search for such events in about 
4000 light curves of microlensing candidates detected by the Optical Gravitational 
Lensing Experiment (OGLE) towards the Galactic Bulge from 1992 to 2007. The 
search reveals a total of 19 repeating candidates, with 6 clearly due to a wide binary 
lens. As a by-product we find that 64 events (~ 2% of the total OGLE-III sample) 
have been miss-classified as microlensing; these miss-classified events are mostly nova 
or other types of eruptive stars. The number and importance of repeating events will 
increase considerably when the next-generation wide-field microlensing experiments 
become fully operational in the future. 

Key words: gravitational lensing - Galaxy: bulge - binaries: general - stars: plane- 
tary systems 



1 INTRODUCTION 

Microlensing events are usually assumed to be single-peaked 
and constant in the baseline before and after the event. Any 
additional brightening episode is usually interpreted as an 
indication of some type of variability other than microlens- 
ing. It is a conservative assumption adopted at the begin- 
ning of microlensing surveys in the early 1990s. Even now 
it is still imposed in microlensing searches in very crowded 
fields heavily contaminated with eruptive variable stars, e.g. 
towards Magellanic Clouds or M31. 

However, microlensing is now well established and many 
events are routinely discovered in real-time. It is well-known 
that in principle a small fraction of microlensing events can 
repeat. By repetition we mean a second brightening episode 
well after the first one such that the star has (largely) re- 
turned to its baseline. A variety of scenarios can cause rep- 
etitions, for example, when the lensed source is a binary 
wide enough for the lens to magnify both components sep- 
arately, one after another. The binarity of source does not 
have to be physical, as the lens can sequentially magnify two 
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unrelated stars blended in one seeing disk. Another mecha- 
nism that can produce two independent peaks in the light 
curve is a lens consisting of two gravitationally bound com- 
ponents with large separation magnifying a single source. 
The lensing of a physical binary by another physical binary 
can result in diverse light curve shapes, which in some cases 
may mimic repeating events. This rather complicated sce- 
nario is not pursued in this analysis of limited number of 
candidate events. One can even imagine two independent 
masses acting as lenses. This is, however, expected to occur 
with a much smaller probability since the microlensing op- 
tical depth for a single event is already low (~ 10 -6 ), and 
thus we do not study such a scenario in this paper. 

iGriest fc Hul (|l992h investigated the effects of the bi- 
narity of sources on the properties of light curves produced 
by a single lens. They concluded that, depending on the 
adopted population of the binary population and the lens 
mass function, the cases with double peaks and asymmetric 
light curves (easily distinguishable from single source events) 
comprise 1% to 3% of all events. 

|Pi Stefano fc Mad <[l996) studied the possibility of re- 
peating microlensing due to wide binary lenses and predicted 
that 0.5-2 per cent of observed microlensing curves should 
exhibit apparent repetitions. The fraction is lower than 
that predicted for "ordinary" binary events with smaller 
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separations (|Mao fc Paczvriskilll99ll '). To date, hundreds of 
ordinary binary microlensing events have been discovered 
llAlcock et all I2OO0I: ICollingel 120041 : Ijaroszvriski et all I2004L 
120061 ; ISkowron et alj|2007t ). but only one "repeating" binary 
event wa s found (OGLE-2q 3-BL G-291) and was studied in 
detail by Ijaroszvriski et al.l l|2005h . 

Thanks to microlensing surveys such a s the Optical 
Gravitational Lensing Experiment (OGLE, UdalskJ 2003 



and M icrolensing Observations in Astrophysics fMOA. lYock 
1998), several thousand microlensing events have been de- 
tected in real-tim^E Among these, we would expect tens of 
events caused by wide binaries. 

Detection and studies of these events are interesting for 
two important reasons. First, these events provide an in- 
dependent way of studying the statistical properties of the 
binary star population, including the distribution of mass 
ratios and separations. Since a large fraction of lens popu- 
lation is expected to be composed of low-mass (~ 0.3Mq) 
main sequence stars, the detection of the binary compan- 
ions around these faint stars would be di fficult using spec- 
troscopic methods fe.g. iFisher et al1l2005l ). because of their 
typically long periods, small accelerations and slow radial 
velocity variations. The novel aspect of the microlensing 
method is that we can infer the mass ratios in the time 
domain simply from the ratio of the timescales in the two 
brightening episodes (see jj5l) rather than using spectroscopy 
(|Di Stefano fc Mag 1 199a ). Second, repeating microlensing 
event s can in principle be used to detect planetary compan- 
ions (Di Stefa no fc Scalz o 1999) ; this would be an extension 
of the methods used currently by the survey and followup 
teams for succes s ful detections of extrasolar planets (e.; 



Bond et al 



Gould et al 



2004 lUdalski et al 1 120051 : iBeaulieu et ail 12006 



2006; lGaudi et alj |2008) 



In this work we present the results of a systematic search 
for repeating events in the available data on all microlensing 
candidates detected by the OGLE team. The paper is organ- 
ised as follows. Sections [2] and [3] describe observational data 
and our search procedure. In Section [4] we focus on different 
microlensing models and algorithms to fit the light curves 
of repeating candidates. The major results of our search are 
presented in Section [5] and we summarise and further dis- 
cuss our results in Section [6] 



2 OBSERVATIONAL DATA 

In this study we use data acquired by the first, sec ond and 
third phases of OGLE. The readers are referred to lUdalskil 
|2003l ) for more technical details about these phases of the 
project. We gather all the microlensing events detected by 
the OGLE Early Warning Sys tem (EWS), by OGLE-II 
(1998 -2000. lUdalski et al.lil994l ) and by OGLE-III (2001- 
2007. IUdalskill2003lL We a lso in clude events found indepen- 
dent lyTylm|niareFanj2TOl|Un the OGLE-II 1997-1999 



data and by Wvrzvkows ki ( 2005) and IWvrzvkowski et al.l 



(|2006h in the OGLE-III 2001-2005 data. In total there are 
4120 unique microlensing candidate events. The systematic 
checking of a possible later (or an overlooked earlier) episode 
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of variability of sources belonging to EWS catalogs is not 
routinely done, and this task is undertaken here. Notice that 
events with closely separated peaks may be removed from 
the candidate list, and thus we focus on "repeating" events 
where the star has (largely) returned to its baseline between 
the episodes since such events will not be rejected by the 
EWS. 

For each event, we check if any additional observa- 
tional data are available in the OGLE-I (1992-1995), OGLE- 
II (1996-2000) and OGLE-III (2001-2007) databases. The 
maximum available time span is therefore 15 years. 152 
microlensing events were observed continuously during all 
three phases of the OGLE project (15 years). About 1200 
events were observed continuously for 10 years and about 
2300 for 5 years. 



3 SEARCHING FOR REPEATING EVENTS 

For each of the 4120 light curves, we search for the presence 
of two or more magnification episodes in the whole time span 
by visual inspection and with a semi-automated algorithm. 
We discuss these two approaches in turn. 

3.1 Visual inspection 

The main criterion used in this subjective analysis is the 
presence of a significant bump before or after the main 
brightening episode. We look for events, where the brighten- 
ings are separated by a return of the light curve to the base- 
line. Our systematic inspection of all the 4120 light curves 
has revealed 13 candidates of repeating events. 

In addition, a small number of objects have been found 
among the microlensing candidates being other types of 
brightenings. They are present only in among events listed 
by the EWS because of the nature of the EWS detection 
algorithm, which detects significant and continuous bright- 
ening, but does not perform any microlensing model fitting. 
The main contamination is from nova-like outbursts and 
from other variables, mostly eruptive stars. The long time 
span allows a better classification of the events, e.g. by iden- 
tifying later outbursts or another period in oscillations. In 
addition most of these stars exhibit some asymmetry in their 
brightening episodes. 

The largest and most uniform sample is from the EWS 
in the OGLE-III phase (until September 2007). Out of 3159 
microlensing candidates in this sample, 64 (~ 2%) turn out 
to be intrinsically variable stars; 24 of these show a be- 
haviour similar to dwarf novae stars with multiple, short- 
time outbursts. Furthermore, 52 events (~ 1.6%) have du- 
plicate entries since they were "discovered" twice in overlap- 
ping adjacent fields. 

3.2 Semi-automated search algorithm 

The main goal for developing an automated algorithm is to 
identify repeating events in which the second magnification 
episode is considerably smaller than the first one, as such 
configuration would be easier to miss visually. Because of the 
contamination by variable stars and caustic crossing events 
we were only able to construct a semi-automated algorithm, 
which still required some human supervision. 
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In the first step we detect t h e mai n microlensing episode 
and fit it with the IPaczvriskH (|l986l ) model. The portion 
of this model light curve which is above 3cr level from the 
baseline defines the duration of the brightening episode, and 
the corresponding data points are removed from the light 
curve. The fit of the Paczynski model does not have to be 
perfect, since it only serves to mark the brightening episode, 
and our experience shows that the same algorithm can also 
be successfully applied to binary lensing events. 

A constant light curve model and the Paczynski model 
are then fit to the remaining data, yielding two x 2 ' s: 
xLd.const, and xLd.Pac- K a repeating event is present, 
X2nd,Pac i s expected to be substantially smaller than 

X2nd, const ■ 

To ensure the second fit is reliable, we calculate the 
number of data points in the second magnification episode, 
ri2nd, and reject all events with «2nd < 3 or ri2nd > N/2, 
where N is the total number of data points. We then con- 
struct the following statistic: 



IX2nd,P ac ~ X2nd .const I 
X2nd,Pac 



(i) 



and choose events with s > 0.2. This ad hoc criterion typ- 
ically corresponds to a fit improvement in Ax 2 of several 
tens, which is high enough to avoid false identifications (cf. 
§5.1.2). On the other hand all 13 events found by visual 
inspection are recovered. 

A total of 193 events have passed the described crite- 
ria. All of them have been examined visually giving 6 new 
candidates, and bringing the total number of candidates for 
repeating events to 19 (see Table [T}. 6 new candidates were 
overlooked during the visual inspection due to their small 
amplitude of the secondary peak. 

3 out of 19 of our candidates were found in pre- 
vious studies: the connection between two EWS events 
OG LE-1999-BUL- 4 2 and OGLE-2003-BLG-220 was found 
by iKlimentowskl ||2005|); OG LE174828.55-221639.9 was 
found by IWvrzvkowskiMgOOSl) and OGLE - 2003- BLG-291 
was described in detail by Jaroszvriski et af] (|2005l ). 



4 MODELING 

All 19 events were fitted with three different models: binary 
source, binary lens and approximate wide binary lens mod- 
els. 

We model a binary lens following ISkowron etail l|2007h . 
In total there are seven parameters. The two point lenses are 
described by the mass ratio (q) and separation (d) in units 
of the Einstein radius (te)- The microlensing geometry is 
described by the impact parameter (6), the angle between 
the source trajectory and the projected binary axis (j3 in de- 
grees), Einstein-radius crossing time (i,E in days), the time of 
the closest approach to the center of mass of the binary (to), 
and the fraction of light contributed to the total blended 
flux by the lensed source, / (/ = 1 indicates no blending). 
The event baseline magnitude (Jo) is measured separately. 
A point-like source is assumed in the fitting. In some cases 
in order to better fit the observational data we further in- 
troduce the parallax motion of the Earth (with a parallax 
scale 7te, defined as la.u./fE, where te is the radius of the 
Einstein ring projected into the observer's plane) and/or the 



rotation of the binary lens assuming circular face on orbits 
(with an angular velocity of /3 in units of deg yr _1 ). 

The first stage of searching for the best-fit models was, 
however, done on a grid more appropriate for uncovering 
repeating events. The grid covers a wide range for 6 param- 
eters: the mass ratio, binary separation, minimum distances 
to the source trajectory from the first and second mass and 
their corresponding times of minimum approach. The start- 
ing values for the approach times are estimated from the 
visual inspection of the light curve. For fixed distances from 
the two masses, two source trajectories are possible depend- 
ing on whether the source trajectory intercepts the binary 
lens axis or not. Each parameter is sampled with about 
15-20 intervals. These initial searching parameters are then 
transformed to the "standard" model parameters mentioned 
above and the \ 2 ls calculated. A few hundred models with 
low x 2 's from the grid are taken as initial guess parameters 
and fed i nto a minimisatio n procedure based on the Powell's 
method (|Press et al.lll992T ). 

We also fit a simple static binary source model to each 
event. In this model, the light curve is a sum of two stan- 
dard single microlensing events with two impact parameters 
(b\ and 62), two parameters for the times of maximum mag- 
nification (toi and £02), fractions of light contributed to the 
total light by the two sources (fi, and /2, /1 +/2 $S 1), base- 
line magnitude (7o) and Einstein radius crossing time (tfi). 
To ensure the resulting models are comparable with the bi- 
nary lens models a similar minimisation strategy is used 
with the same grid sizes and optimisation routines as those 
for the binary lens model. Notice that the binary source 
model was fit only to non-caustic crossing events since the 
binary source models cannot reproduce the sharp gradient 
features in caustic-crossing events. In addition, we apply 
an approximate wide binar y model following the concept 
of iDi Stefano fc Mad | 19961 ). In this model, the binary lens 
acts as two independent single lenses. If the two individual 
magnifications are given as [11 and /i2, then the resulting 
magnification is approximated as fj, w fix + fj,2 — 1- The fit- 
ting is done on a grid similar to that for the full binary 
lens model. In practice it is important to fit both bumps 
simultaneously to ensure the same constraints on the blend- 
ing (and other) parameters. Due to the degen eracies in the 
blending model (|Wozniak fc Paczynski 1997), fitting each 
peak independently gives uncertain estimation of fluxes and 
timescales. 

Notice that the two (smooth) peaks produced by a wide 
binary lens correspond to the source approaches to the di- 
amond shaped caustics, not to the lens components. Since 
both caustics lie between the masses, fitting the simplified 
model systematically under-estimates the binary separation 
and may give incorrect trajectory direction. However, since 
we are primarily interested in the mass ratio of the binary 
which are given by the rati o of the squares of the tim escales 
of two magnification peaks (|Di Stefano fc Maolll996h . this is 
not a significant limitation (see 33 and Fig. [2| . 



5 RESULTS 

Table Q] lists all 19 candidates for repeating events. It shows 
also the mass ratios in the full and approximate wide bi- 
nary lens models and the light ratio from the binary source 
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Table 1. Candidates for repeating microlensing events. 



Event OGLE- 



{h/h)hs 



<3bl <?approx Type 



1999-BUL-42 




model not found 






cc bl? 


1999-BUL-45 





189 


0.203 





270 


bl/bs 


2000-BUL-42 





260 


0.339 





400 


bl 


2002-BLG-018 





280 


0.395 





443 


bs/bl 


2002-BLG-045 





045 


0.008 





011 


bl 


2002-BLG-128 






0.611 






cc bl 


2003-BLG-063 





539 


0.203 





428 


bl 


2003-BLG-067 t 





817 


0.788 





844 


bs/bl 


2003-BLG-126 





292 


0.604 





918 


bl 


2003-BLG-291 






0.617 






cc bl 


2003-BLG-297 





075 


0.147 





121 


bl/bs 


2004-BLG-075 





707 


0.587 





664 


bl/bs 


2004-BLG-328 





150 


0.056 





054 


bs 


2004-BLG-440 





927 


0.622 





906 


bl 


2004-BLG-591 





548 


0.507 





431 


bl 


2006-BLG-038 






0.569 






cc bl 


2006-BLG-460 





388 


0.267 





285 


bl 


175257.97-300626.3 






0.195 






cc bl 


174828.55-221639.9 






0.158 






cc bl 



Light ratio from binary source model and mass ratios from full 
(%i) and approximate (q a pprox) binary lens models are shown. 
The last column indicates whether the binary lens ("bl") or bi- 
nary source ("bs") model is better, "bl/bs" stands for comparable 
binary lens and binary source models. Events with clear caustic 
crossing features ( "cc bl" ) admit only the full binary lens model, 
f Modelling for this event indicates it did not fully return to the 
baseline between the peaks (see the light curve in the Appendix) . 
Nevertheless, we include this event because of the lack of sufficient 
data points around the "saddle point" between the peaks. 



model. The last column in the Table indicates the nature 
of each event as concluded from a comparison of the best 
X 2 values of the binary lens and binary source fits. The two 
models are regarded as comparable when their \ 2 values 
differ by less than 10. In our terminology, "bl" stands for bi- 
nary lens, "bs" for binary source and "cc bl" indicates that 
the caustic crossing features are clearly visible in the light 
curve for which, therefore, only the full binary lens model 
is fit. The parameters of the best-fit binary lens and binary 
source models are listed respectively in Tables [2] and [3] The 
light curves of all the events are shown in the Appendix. 

The smooth (non-caustic-crossing) light curves with two 
peaks quite often have concurrent, similar quality fits based 
on binary source and binary lens models. This degeneracy 
is also known for events wit h much shorter duration , when 
the peaks partially overlap jjaroszvriski et all |2004| . I2006J ; 



ISkowron et al.ll2007l ; ICollinge]|2004 ). Since binary lenses offer 
a much larger variety of possible light curve shapes, their fits 
are usually formally better. 

Clear caustic crossing features can be seen in the light 
curves of about 30 per cent of our candidates. This frac- 
tion is an order of magnitude higher than that in the whole 
sample of microlensing events - in the uniform sample of 
3159 candidates from the EWS from OGLE-III, 73 events 
(~ 2.4%) show clear caustic crossing features in their light 
curves. Our simulations (section l5.1|l show that assuming all 
our candidate events were caused by wide binary lenses we 
should expect that approximately 10-30 per cent should ex- 
hibit caustic crossing features (the fraction increases as the 
mass ratio decreases). The high fraction of caustic crossing 
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Figure 1. Time between two peaks in repeating events versus 
the Einstein radius crossing time, tg. 



events in our sample of repeating events suggests that the 
binary lens scenario is favoured for the majority of our can- 
didates. 

The number of binary source events in our sample ap- 
pears to be too low when compared w i th predict i ons o f 
iGriest fc Hul (|l992T l. Irian fc Jeond rtl998l ). iDominikl l|l99ct ). 
Even if all ambiguous cases were classified as binary source 
events, we would get probability of only 0.15% for their oc- 
currence, much lower than the percentage from these stud- 
ies. Our sample, however, includes only the events with well 
separated peaks, while the other authors include close binary 
sources, which have the highest probability of producing an 
event distinguishable from ordinary microlensing light curve. 

One practical issue one needs to know for future search 
for repeating events is how long an observer has to wait 
for the secondary episode to occur. This is illustrated in 
Figure [TJ showing the time betw een the two peaks vs. the 
event timescale. As predicted in |Pi Stefano fc Maol (|l996l ) 
this time is of the order of a few (from 2 to 6) Einstein-radius 
crossing times. For the repeating events described here this 
is between 32 to 472 days with a median of 142 days. 

Since the approximate binary model gives a quick way 
to estimate the binary mass ratios (see it is important 
to check its accuracy. Figure [2] compares the mass ratios in 
the full and approximate binary lens models where both fits 
are available. A strong correlation is clearly visible (with an 
RMS difference of about 36%), indicating that the simpli- 
fied model can indeed be used to estimate the mass ratio in 
repeating events due to wide binary lenses. 



5.1 Distribution of binary lens mass ratios 

The histogram in the bottom panel of Fig. [2] shows the dis- 
tribution of the mass ratio in the wide binary lens models. 
It shows a peak around q ~ 0.5 and a decrease for small 
values of q. The decrease is, however, likely due to a larger 
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Table 2. Parameters of the best-fit binary lens models. 



Event 


X 2 /DOF 




1 


d 







b 


to 


t 


E 




/ 




Io 


1999-BUL-42 








2003-BLG-220, 


model not found, see text 












1999-BUL-45 


645.4/405 


0. 


203 


4.104 


18. 


00 


0.47 


1401.6 


30 


.2 


1 


.00 


17 


.75 


onnn t>ttt a o 

2000-BUL-42 


8852.0/697 


0. 


339 


5.308 


22. 


10 


1.09 


2096.7 


99 


.0 





.37 


13. 


.58 


2002-BLG-018 


142.9/111 


0. 


395 


6.572 


-1. 


.30 


0.44 


2510.3 


34 


.6 





.91 


18. 


.06 


2002-BLG-U45 


3222.7/633 


0. 


008 


3.958 


183. 


49 


-0.22 


2360.8 


26 


.4 





.97 


18. 


.76 


2002-BLG-128 


2218.6/810 


0. 


611 


1.908 


173. 


84 


0.22 


2463.2 


59 


.5 





.12 


17 


.74 


onnQ dt n^Q 
ZUUo--t>-LLr-UOo 


line; /t /Qi7 
11U0.4/O1 f 


u. 


2Uo 


3.154 


25. 


.76 


1.39 


2869.3 


54 


.0 


U 


. 1 O 




1 Q 


2003-BLG-067 


1236.4/323 


0. 


.788 


3.636 


0. 


.60 


0.47 


2942.4 


88. 


.6 





.62 


16. 


.45 


2003-BLG-f26 


1027.9/295 


0. 


604 


2.808 


-34. 


80 


-0.92 


2804.2 


22 


.8 


1 


.00 


15 


.66 


2003-BLG-29ft 


912.8/251 


0. 


617 


3.041 


184. 


70 


0.50 


2925.8 


13 


.5 





.38 


17 


.45 


2003-BLG-297 


984.6/297 


0. 


147 


2.912 


179. 


33 


0.51 


2963.3 


75 


.5 


1 


.00 


17. 


.31 


2004-BLG-075 


555.5/203 


0. 


.587 


4.456 


-11. 


.09 


-0.19 


3097.8 


9 


.7 


1 


.00 


18. 


.82 


2004-BLG-328 


2782.3/545 


0. 


056 


3.399 


175. 


70 


0.20 


3181.5 


25 


.3 





.20 


18. 


.45 


2004-BLG-440 


2220.0/527 


0. 


622 


4.499 


158. 


66 


0.65 


3183.1 


8 


.6 





.36 


16 


.34 


2004-BLG-59f 


2478.3/672 


0. 


.507 


2.874 


4. 


.69 


0.26 


3426.9 


81 


.1 





.29 


18. 


.18 


2006-BLG-038 


3556.7/744 


0. 


.569 


3.138 


206. 


.47 


-0.68 


3809.4 


12 


.7 


1 


.00 


16. 


.44 


2006-BLG-460 


2418.9/862 


0. 


.267 


3.354 


178. 


85 


0.11 


3982.7 


20 


.0 





.69 


19. 


.10 


i75257.97-300626.3t 


4845.1/932 


0. 


195 


3.268 


-171. 


.58 


-0.37 


2547.4 


102 


.9 





.35 


18. 


.19 


174828. 55-221639.9§ 


406.2/192 


0. 


.158 


2.431 


172. 


.48 


-0.06 


2827.7 


155 


.2 





.31 


18. 


.93 



The columns are respectively, the name of the OGLE event, x 2 / rmm t , er of degree of freedom (DOF), mass 
ratio q, binary separation d (in units of the Einstein radius), direction of the source trajectory with respect 
to the binary axis fi (in degrees), impact parameter b, time of the closest approach to the center of mass fo 
(Julian date shifted by 2450000), Einstein radius crossing time tg (in days), blending parameter / = F B /Fo 
and I-band brightness of the baseline io- 

f The model for the event 2003-BLG-297 is taken from Jaroszynski et al. (2005) which includes binary 
rotation. | model with parallax motion, tte = 0.279. § model with binary axis rotation, = 0.032° /day. 



Table 3. Parameters of the best-fit binary source models. 



Event 


X 2 /DOF 




61 




b 2 


toi 


t02 


t 


E 




fl 




h 




Io 


1999-BUL-42 


































1999-BUL-45 


646.4/405 





.450 


0. 


.546 


1310.64 
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The columns show the name of the OGLE event, x 2 A ne number of DOF, impact parameters 61 and b 2 for 
two source stars, times of closest approach to both components toi and U12, Einstein radius crossing time tE, 
blending parameters /1 = F s i/(F s i + F s 2 + Ff,) and f 2 = F s 2/(F 3 i + F s2 + Fj,), and I-band brightness of 
the baseline Io- F s i and F s 2 are the fluxes from the two sources, and Ft, is the flux from any other unrelated 
star(s) within the seeing disc (blend). 



incompleteness at small q in the observations. Thus to un- 5.1.1 Simulations of detection efficiency for wide binary 
cover the true distribution of mass ratios, we have to correct lenses 
the histogram in Fig. [5] for the detection efficiency. 



A rigorous detection efficiency analysis at the pixel level 
involves many factors and is beyond the scope of the paper. 
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0.01 0.1 1 

q 

Figure 2. Correlation between the mass ratios obtained in the 
full (q) and approximate (g a pprox) binary lens models (upper 
panel) and the observed histogram of mass ratios derived in the 
full binary lens model (lower panel). 



However, since our sample is small and Poisson errors are 
large (see Fig. [3j , a first-order approach at the catalog level 
is sufficient to account for the primary selection effects for 
wide binaries with different mass ratios. 

We perform simulations of repeating events due to wide- 
separation lenses with the mass ratio ranging from 10 -3 to 1 
at 25 uniform, logarithmically-spaced intervals. Depending 
on the mass ratio, about 10 4 — 10 5 binary lensing events 
are generated. An automated algorithm for identifying the 
fraction of repeating events observed in the sample is then 
used (see £15.1. 2[) . 

To create a synthetic microlensing event, we first gener- 
ate a physical model of a binary lens and source trajectory 
and then based on the physical model, a light curve with 
sampling rates similar to those of OGLE is simulated. To 
generate a binary microlensing model the following steps 
are performed: 

(i) For a given mass ratio, a binary separation is drawn 
from the range 1 to 36 Einstein radii using a uniform dis- 
tribution in logarithm. Our simulations show that the de- 
tection efficiency falls to zero outside this range. For lenses 
of masses M £ [0.1, 1] M© placed in the Galactic disk and 
sources in the bulge, the above range corresponds roughly 
to physical separations of 1 - 100 AU. 

(ii) To ensure at least one magnification episode, we set 
the source trajectory to pass the heavier lens at an impact 
parameter from to 1 Einstein radius. 



(iii) The trajectory angle is set randomly from to 360 
degrees. 

(iv) The time of the closest approach to the heavier mass 
is chosen from a uniform distribution to lie within the whole 
observing time. 

(v) The Einstein-radius crossing time is drawn from an 
approximate distribution that matches that of the observed 
OGLE lenses: a log-normal distribution centered on tE = 26 
days with a standard deviation of 0.3. 

After the physical model is chosen using the procedure de- 
scribed above, the magnification as a function of time can 
be easily found. However, to mimic an observed light curve, 
several further steps are necessary. 

(i) The baseline magnitude is chosen randomly from the 
luminosity function of one of the OGLE fields (BLG104.6). 

(ii) The value of the blending parameter is chosen at ran- 
dom bet ween and unity. The small sample of binary lens 
models (|Jaroszvriski et al]|2004l ) supports such a choice of 
the blending parameter distribution. The values of the base- 
line flux and blending parameter combined with the theo- 
retical lens model yield a perfect light curve without errors. 

(iii) To account for observational gaps in the light curve, 
the epochs of observations are chosen so as to match ei- 
ther BUL J SC34.I.36546 in OGLE-II or BLG104.6.I.7723 in 
OGLE-III. They have been chosen as a reference for epochs 
and observational errors (to be used in eq. [2| as their light 
curves span from 1998 to 2007 during which all our repeating 
events candidates have been found. 

(iv) Each measurement has an observational error AI 
from the rescaled value s for a reference star A/ re f using the 
empirical formula from I WvrzvkowskH l|2005h : 

AI = AJ rcf 10°' 33875(/ - 7 " f) (2) 

where I is the model magnitude, I Te f is the magnitude and 
A7 re f is the error bar for the reference star at the epoch. 

(v) Gaussian errors are added to the light curve. In this 
step, we take into account the fact that OGLE errors are 
slightly underestimated: for every simulated data point, a 
Gaussian standard deviation is derived based on its error 
bar value AI in eq. ((2]) using 

a = aAi.SSA/) 2 + 0.0052 2 . (3) 

This rescaling has b een o btained following the method of 
IWvrzvkowski et al.l (|2008l ) by comparing the rms of a set 
of constant stars to their mean error bars returned by the 
photometry pipeline. 

5.1.2 Efficiency- corrected statistics of binary mass ratios 

To calculate the detection efficiency of repeating events for a 
given mass ratio, all the light curves of synthetic events must 
be classified. To accomplish this, we fit a standard Paczyhski 
model to all the light curves and compare the resulting x 2 
with a constant line fit. We have also applied the same pro- 
cedure to check if a Paczyhski model fit to a constant lu- 
minosity light curve (with Gaussian errors) can improve the 
fit. Our experiment shows that the typical improvement for 
~ 10 3 d.o.f. is Ax 2 » 15 and exceeds 55 in - 1% of cases. 
When investigating the synthetic light curves, we require 
a A\ 2 > 55 improvement of the Paczyhski model over the 
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Figure 3. Relative detection efficiency, defined as the ratio of 
the number of detected "repeating" events and the total number 
of detected microlensing events, as a function of the mass ratio 
q. Error-bars denote Poisson noise based on population of simu- 
lated sample. The solid line is an empirical fit of N rep /N cvcnt fa 
0.021g - 687 . 



constant line model to treat the event as a microlensing can- 
didate. (Notice that while each curve is simulated using mi- 
crolensing, the gaps between the seasons and/or sampling 
rate may prevent it from being "discovered".) 

We then use our automated algorithm (see i]3.2|l to iden- 
tify repeating events in our sample. The relative detection 
efficiency is then defined as the ratio of repeating events 
(N Icp ) found in the sample of all identified microlensing 
events ( Advents)- Figure [3] shows the efficiency with respect 
to mass ratio (q). It can be approximated by the formula 

iVrcp/iVcvcnt ~0.021 (/ °- 687 . 

Fig. [3] shows likely underlying distribution of binary 
mass ratios obtained by convolution of observed mass ratios 
(Fig. [2} and their detection efficiency (Fig. [3}. For mass 
ratios between 0.07 < q < 1, the distribution is consistent 
with a uniform distribution in the logarithm of mass ra- 
tio with Poisso n errors, in agreement with earlier results of 
I Trimble! ^199($ ). Although we have one event (OGLE-2002- 
BLG-045) for the smallest mass ratio bin, because of the 
low detection efficiency for such events, the inferred relative 
number is quite high, but the error bar is large. Neverthe- 
less, it may indicate the presence of a different population 
of binary (planetary) systems with extreme mass ratios. 

Assuming a uniform distribution of binary lenses in 
log q and log d, we can calculate the averaged efficiency of 
the detection for q € [0.1,1] and d G [l,36]r_e, obtaining 
< A^cp /Agents >= 0.0105. Notice that we have excluded 
the brown dwarf and planetary binary companions (with 
mass ratios between q 6 [0.001,0.1] since there is only one 
such planetary candidate. If all the events were caused by 
binary lenses with parameters from this range, one would ex- 
pect 4120 x 0.0105 ~ 43 detections of repeating events. Since 



i i i | 1 1 — i — i — i i i i | r 

2500 ~ 7 T 

2000 r 

S-H 



,£1500 - 




0.01 0.1 1 

q 

Figure 4. Distribution of binary mass ratios after correcting for 
detection efficiency. Error bars are from the Poisson noise. The 
left-most bin at q = 0.01 contains only one event. 



twelve repeating events are unambiguously classified as bi- 
nary lens cases, this suggests that at least 28% of all lenses 
are binary systems within the specified range in the mass 
ratio and separation. Binary stars roughly follows a uniform 
distribut ion per decade of separation bet ween 10 cm to 
10 17 cm l| AbtJl 19831 ; iMao fc Paczvriskilll99ll) . thus we expect 
a fraction of log(36/l)/6 = 26% within d € [l,36]r_B, in good 
agreement with the 28% we estimated above if all the lenses 
are in binary systems. However, earlier studies indicate that 
about 50% of solar-type stars are in bi nary stars, while for M 
stars this fraction may be even lower uennov fc M ayor 
Il99ll ; iFisher fc Marcel 19921 ; iReid fc Gizislll997l ; lLadalbOOrj 
and references therein), thus we would expect a factor of 
~ 2 smaller number of repeating events due to binary lenses. 
There are several possiblities to explain the "discrepancy". 

1) It may simply be due to small number Poisson fluctuation. 

2) The binary fraction is under-estimated in spectroscopic 
surveys. 3) The distributions of the separation and mass ra- 
tio may differ from being uniform per decade as we assumed 
(our data are consistent with this assumption within the 
small number statistics.) 4) Some of the "repeating" events 
are actually not due to binary lenses. In the future, with 
a much larger number of repeating events and denser time 
samplings (e.g., with OGLE-IV), we will be in a better po- 
sition to test these different possibilities. 

5.2 Comments on two individual events 

In this section, we study two particularly interesting "re- 
peating" microlensing events. OGLE-1999-BUL-42/OGLE- 
2003-BLG-220, and OGLE-2002-BLG-045 in more detail. 

5.2.1 OGLE-1999-BUL-42/OGLE-2003-BLG-220 

For this event, we were not able to find any satisfactory 
model. This interesting case requires closer scrutiny. The 
light curve of the event in 1999 resembles a caustic cross- 
ing binary lens, but because of the lack of sharp features it 
can also be explained by a binary source model. However, a 
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static binary source model turns out to be insufficient. The 
secondary bump from 2003 shows nearly no binarity features 
except one outlying data point and a slight asymmetry. 

The timescales of both magnification episodes are of the 
order of 20 days and the separation between them is very 
long, about 1400 days. This indicates it is unlikely that both 
bumps were produced by a very wide binary lens, because 
such lenses have tiny caustics and as such they cannot ex- 
plain the features in the light curve from 1999. One possible 
scenario for the repeatability is that the source was lensed 
by two unrelated lenses or two sources were lensed by one 
binary lens; the event may also be due to a variable star 
rather than microlensing. 

We have attempted to gain some additional insights 
using information available from the OGLE photome- 
try pipeline u sing the Difference Image Analysis (DIA) 
l|Wozniakl l2000) . In this method all variable sources are de- 
tected on the subtracted image and for each variable the 
centroid of the subtraction residuals indicates the real posi- 
tion of the variable star. Because of blending this position 
may not be aligned with the position of the baseline star 
on the template image. In the OGLE database (A. Udal- 
ski, private communication) we checked the positions of the 
residuals with respect to the template positions around the 
1999 and 2003 peaks. Because of sparse sampling there are 
only a few high signal-to-noise measurements in the 2003 
peak, however we can still confirm that the position of the 
lensed source is consistent with the centre of the blend. This 
does not necessarily mean there was no blending, but sim- 
ply that the source happens to be located near the centre 
of the blended object. On the other hand, the 1999 peak 
had a much better coverage and the detected astrometric 
signal is firm, showing a clear displacement of the residuals 
in respect to the blended position by about 150-200 mas. 
This result indicates that it is unlikely the 1999 and 2003 
events were caused by the same lens, unless it moves very 
fast (~ 40mas/yr) and/or is a very close object. The re- 
maining (microlensing) hypothesis is that it is a pair of com- 
pletely unrelated events occurring on two different sources. 
Further studies, including high resolution imaging (e.g. with 
the Hubble Space Telescope) , would either confirm or reject 
this scenario. 



5.2.2 OGLE-2002-BLG-045 

The second brightening of this event was much shorter than 
the first one and the binary lens model gives a low mass 
ratio of q — 0.008 (the approximate binary lens model gives 
q = 0.011). For a typical lens mass of about O.3M0, this 
would imply a planet mass of a few Jupiter masses. Un- 
fortunately, sparse sampling of both peaks prevent us from 
concluding convincingly the planetary nature of the lens. 
Nevertheless this demonstrates repeating events may be a 
potentially exci ting channel of detecting p lanets in repeat- 
ing events fe.g-. lDi Stefano fc Scarzolll999t) . 

5.3 Contaminants to the OGLE Early Warning 
System 

As a by-product of our search, we identified a number of 
non- microlensing events in the database. In total we found 



64 mis-classified events out of 3159 events (about 2%) in the 
OGLE-III microlensing candidates identified by the EWS. 
This low rate is re-assuring as it indicates that EWS, par- 
tially based on human interpretation, is highly reliable in 
discovering true microlensing events. Dwarf novae are the 
main contaminators (24). Thanks to the 15 years of con- 
tinuous monitoring of OGLE, it is possible to detect up to 
30 outbursts for one dwarf nova. However, the majority of 
our stars had 2-4 observed outbursts. It is also possible to 
measure very long pseudo-periods of outbursts - for 5 of our 
dwarf novae candidates the time between outbursts is longer 
than 1000 days. Note that all dwarf novae found with the 
visual inspection are also retrieved by our automated fitting 
algorithm. Thus it is possible to use this approach to iden- 
tify these stars in the entire OGLE photometric database; 
we plan to conduct such a study in the near future. 



6 CONCLUSIONS 

A small fraction (~ 0.5%) of microlensing events do repeat: 
our search revealed 19 repeating candidates in the sample 
of 4120 microlensing events we studied. This gives a rate of 
~ 3 events per year in the OGLE database at the current 
discovery rate of ~ 600 events/yr. Both the high fraction 
of caustic crossing events and comparison of the goodness 
of fit between binary lens and binary source models suggest 
that probably most of these are due to wide binary lenses, 
although the predicted number seems to be somewhat higher 
than expected (see H5.1.2|) . 

We find that it is possible to estimate the mass ratios for 
repeating events in a straightforward manner. With a grow- 
ing number of microlensing events observed every year these 
events could provide a valuable sample to study the stellar 
binary populations in the Galaxy. This method operates in 
the time domain, different from the usual spectr oscopic stud- 
ies of binaries in the solar neighbourhood (e.g.. iFisher et all 
2005). Accounting for the detection efficiency, the distribu- 
tion of mass ratios of binaries appears to be consistent with 
a uniform logarithmic distribution, i n agreement w ith that 
from previous spectroscopic studies |Trimbkjll990t ). 

Our study also illustrates an example of a missed oppor- 
tunity to find a planet (OGLE-2002-BLG-045). The second 
brightening episode was much shorter than the first. Unfor- 
tunately the sampling was not dense enough after the first 
episode returned to the baseline. In the future it would be 
profitable for survey teams and follow-up networks to pay 
more attention to microlensing events even after the main 
magnification peak (the median time between the two peaks 
is about 5 months). However, given the limited observational 
resources, this is difficult to implement in the current mode 
of extrasolar planet discovery with alerts from survey teams 
followed by other teams with intensive observations. Nev- 
ertheless, the number of repeating events will increase con- 
siderably with next-gener ation wide-field m icrolensing sur- 
yeys from the ground (e.g. lGould et alj 2007) and from space 
l|Bennett et al.l2007l ). The ground- based experiments are al- 
ready evolving toward a wide-field network, starting with the 
upgraded MOA-II experiment and the soon-to-be upgraded 
OGLE project (OGLE-IV) . The dense sampling will be par- 
ticularly important for the detection of extrasolar planets on 
wide orbits and will offer a new channel for extrasolar plan- 
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ets discovery (|Di Stefano fc Scalzo|[l999h . For these events 
the mass ratio can be approximately "read" off from the 
light curve using the approximate binary lens model. When 
the next-generation microlensing experiments become fully 
operational, this method may become fruitful. 
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APPENDIX A: LIGHT CURVES OF 19 
REPEATING CANDIDATE EVENTS 
TOGETHER WITH THE BEST-FIT BINARY 
LENS MODEL. 
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